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Abstract 

Spontaneous symmetry breaking (SSB) for Bose-Einstein condensates cannot treat phase off- 
diagonal effects, and thus not explain Bell inequality violations. We describe another situation 
that is beyond a SSB treatment: an experiment where particles from two (possibly macroscopic) 
condensate sources are used for conjugate measurements of the relative phase and populations. 
Off-diagonal phase effects are characterized by a "quantum angle" and observed via "population 
oscillations", signaling quantum interference of macroscopically distinct states (QIMDS). 
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If two or more Bose-Einstein condensates (BEC) merge, they produce an inter 
pattern in their densities, as shown by spectacular experiments with alkah atoms jl]. The 
usual explanation assumes spontaneous symmetry breaking (SSB) of particle number con- 
servation, where each condensate gets a (random) classical phase and a macroscopic wave 
function: 



(V'„,,(r)) = 7n„,,(r) e^^-^W (1) 



where na^pir) and (f)a,i3 are density and phases of condensates a, (3. Alternatively, one can 
use a "phase state" describing two condensates with a relative phase and a fixed total 
number of particles: 

\<P,N) = -^^{ai + e^^aX\0) (2) 

where aj^ and create particles in condensates a and (3, respectively. However, one can also 
consider that two condensates are more naturally described by a double Fock state (DPS), 
a state of definite particle numbers, for which the phase is completely undetermined: 

\No.N,) = -^J== ^^4^^ |0) (3) 

It is found that repeated quantum measurements of the relative phase of two Fock 

states can make a well-defined value emerge spontaneously, but with a random value. For 
example, the probability of finding M particles, out of a total of N, at positions ri. . . tm 
where M ^ is shown to be given by js, ^: 



P(ri, ■ ■ ■ r^) ~ n [1 + cos(k ■ r, + A)] (4) 



r- dX 

i=l 

Positions can be obtained one by one from this distribution; for large enough M the integrand 
peaks sharply at a single value, just as a particular phase is found in the interference 
measurement of Ref . jl| . 



One can ask whether the SSB approach is appropriate js] and whether it gives complete 
information [4]. Indeed we will show that the assumption that the condensates are described 
by Eq. gives a broader range of physical possibilities, which are unavailable when using 
Eq. (|2]). The additional effects involve phase ojf -diagonal terms, which can result in (I) 
violations of local realism, i.e. violations of Bell inequalities, and (II) the occurrence of 
quantum interference between macroscopically distinct states (QIMDS), as discussed by 
Leggett Ql- Neither of these effects is available in the SSB treatment. We have previously 



discussed [4] violations of Bell inequalities with double Fock states. Here we will show that 
the effect II can be detected in interferometer experiments by the observation of "population 
oscillations", first introduced by Dunningham et al 8| within a three-condensate position 
interference analysis. Theses oscillations are more robust than Bell inequality measurements, 
since a few missed particles can be tolerated. 

Leggett 0] considers how one might test for QIMDS by finding coherent superpositions 
involving large numbers of particles ( "Schrodinger cats" ) . One can tell the difference between 
such a pure state and a statistical mixture of the elements of the state only by observing the 
off-diagonal matrix elements between the different wave-function elements. For example, in 
a state of the form \E' = Caipa{^, 2, 3, ■ ■ ■ , A^) + Cbipb{^, 2, 3, ■ ■ ■ , A^) one hopes to see terms 
like {ipa\G\ipb) and its complex conjugate, where G is an appropriate A^-body operator 
connecting the two states. As Leggett j3] says, " ■ ■ ■ what matters is that not one but a 
large number of elementary constituents are behaving quite differently in the two branches." 
Here we discuss an experiment where particles from each of two Bose condensate sources are 
either deviated via a beam splitter to a side collector or proceed to an interferometer. The 
measurements in the interferometer create the two branches, and the detection in the side 
detectors (involving the connecting G operator) allows the observation of the off-diagonal 
matrix elements of the two components. In recent years several experiments have begun to 
make progress toward the goal set by Leggett, by use of large atoms jol, superconductors jlo| . 



magnetic molecules a quantum dot "molecule" 12 1, and photons [13|, and including a 
Bell inequality violation in a Josephson phase qubit 14 1. 

Fig. [U shows the interferometer. The QIMDS state is created from the double Fock 
state by an interference measurement at beam splitter BS, with detectors 1 and 2 giving 
results mi, m2- The path difference between the two sources to BS is represented by angle 6. 
Detectors 3 and 4 record nia and particles respectively; although they seem to measure 
only the source populations, they are actually sensitive to QIMDS, as we will see. 

With a single quantum particle crossing two slits, which act as sources giving rise to 
interference, one can measure either the interference pattern and have access to the relative 
phase of the sources, or from which source the particle comes; they are exclusive measure- 
ments. Here, because condensates provide many particles in the same quantum state, some 
of them can be used for a phase measurement, others for a source measurement. 

The destruction operators 01,02,03 and 04 associated with the output modes of the 
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FIG. 1: Two source condensates, with populations Na and Np, emit particles. Some of them reach 
the central beam splitter BS, followed by detectors 1 and 2 registering mi and m2 counts. The 
other particles are then described by a quantum superposition of macroscopically distinct states 
propagating inside the region shown with a dotted line; they eventually reach counters 3 and 4, 
which register nia and counts respectively. A phase shift 6 = tt/2 occurs in one path. 

interferometer can be written in terms of the source mode operators, and ct/j, by tracing 
back from the detectors to the sources, with a factor 1/a/2 at each beam splitter and a phase 
shift of 7r/2 at each reflection: 



Oi = - {e'^aa + iapj ; 0.2 = - {ie^'^aa + ap 



% t 

The probability amplitude for finding particle numbers {mi, m2, m^} is: 



ma „"^/3 «"*2 

tig Qj^ tfcl (X'2 



^mi\m2lma^.mfs\ 

The double Fock state (DFS) \NaNj3) can be expanded in phase states as 



) (6) 



(7) 



where the phase state having constant total numbers of particles is given by Eq. ([2]). These 
states have the property that, for = fjafla + fi/^a^ {i = 1, 2): 



a, 10, N) = y^(t;,„ + v./se'"^) |0, N - 1) 



(8) 
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so that the state created by the interferometer is: 

|r) = ar^r \NaN^) - 
where M = mi + 1712 and: 



iNf 



n 2n 



N-M) 



(9) 



If we take 9 = % /2 (as we do henceforth) this takes the simple form: 



(10) 



i6/2\M 



COS — sm — 
. 2/ I 2 



(11) 



Fig. Elshows i?(0) = i?(0) X (2ze 



i<f>/2\-M 



which has two peaks at ±0o = =1=2 arctan J 1712/ mi 



1.0 
0.5 



-1.0 \ /-0.5 



0.5 1.0 



-0.5 
-1.0 

FIG. 2: Variations of i?(0) if results mi = 17 and 772-2 = 83 are obtained. The peaks are at 
00 = iO.TSvr (the phase choice 9 = tt/2 gives symmetrical peaks about zero). The relative sign 
of the two peaks is (—1)"^^. For large numbers of particles, the measurement produces a coherent 
superposition of macroscopically distinct states ( "Schrodinger cat"). 

This is not surprising: classically, the ratio of the intensities in the output arms of the 
interferometer determines the absolute value or the phase difference between the two input 
beams, but not its sign. Separating the negative and positive contribution of provides: 



|r) = |^+) + (-i)'"^ 1^-) 



(12) 



Let us begin with a qualitative calculation. We assume that M is large, so that the peaks 
are sharp and: 

|V;±) ~ e^'(^''-*'/2)'^« I ±00, N-M) (13) 

These two wave-function branches are orthogonal for large M for any 0o not too near zero; 
and they are macroscopic as long as — M is large. 



Showing QIMDS requires making a measurement that is sensitive to the interference 
between the two components; this is the role of the side-detectors shown in Fig. [TJ Because 
a^^a^" \±(t)o,N - M) ~ e^''"'^'^'' |0), the probabihty of getting the set {mi, ma, m^, m/j} 
becomes: 

P(mi, m2, m^, mp) ~ 1 + (-l)™^ cos [(m„ - mp) 0o] (14) 

(if Na = Njb), where the cosine terms arises from the sum of the two cross terms 
V'lt a^3^" a]^'^ a'^'^ a^'^ V't)- Now, if one does the interferometer experiment for fixed source 



numbers, say, = Np, and considers only those experiments having the same mi, 1712, then 
the interference between the two elements will show up in a cosine variation of probabilit 
with rria. We call this effect "population oscillations"; it was already discussed in Ref. |6 
for three-condensate experiments. 

These oscillations are beyond SSB since they disappear if one starts from either ([1]) or 
(j2]). With a phase state of phase x for instance, the action of the destruction operators ai^2 
on this state introduces x instead of an integration variable into (|H]); this leads essentially 
to dHj) without the (p integral. No interference effect between two phase peaks occurs and the 
probability is proportional to |-R(x)P- One gets a rria, dependence of the probability that 
is proportional to a simple binomial distribution {N — M)\ / mj.mpl, without any oscillation. 
Actually the angle x plays no role at all in this dependence, which is natural since detectors 3 
and 4 do not see an interference effect between two beams; they just measure the intensities 
of two independent sources after a beam splitter at their output. 

A more accurate calculation is now presented. Operating on Eq. ^ with a™" a™'', and 
forming the probability introduces another angle 0', so that the probability for finding the 
set {mi,m2,ma,mp} takes the form: 



mi!m2!ma!m_g! 2^+*^ J-tt 2tx J-w 271 

(15) 

We note that one phase branch peak occurs for — vr < < and the other for < < vr, 
so that the overlap between different branches occurs for 0' 7^ 0. A change of variables to 
A = (0 - 0')/2 and A = 7r/2 - + (0 + 0')/2 leads to: 

P(mi,m2,m.,m,) = ^^,^^,;;,j^^,2^ ^ ^os [{N^ - - N, + m,) A] - 

X [cos A + cos A] [cos A - cos A] (16) 
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When 6 = Ti/2 the "classical phase angle" A is half the sum of and 0'. The expresssion 
also contains another angle, A, which we call the "quantum angle" - in Ref. {4| it appeared 
as a consequence of a conservation rule, but here we introduce it to characterize quantum 
interference effects between different values of the phase. 

To examine the behavior of the probability, we plot the quantity: 



F(A, A) = [cos A + cos A]""' [cos A - cos A] 



n M-mi 



(17) 



We take A''^ = = M = 100 in our examples here. F{A, A) has multiple peaks as shown 
in Fig. |3]for mi = 17. The extrema are easily shown to occur at: 
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FIG. 3: Plot of F{A, A) as a function of A and A for mi = 17 and m2 = 83. The peaks along A = 
and ibvr correspond to phase diagonal matrix elements, while the negative depressions, having 
A 7^ 0, correspond to off-diagonal matrix elements between two macroscopic phases (QIMDS). If 
mi and ni2 are even, the negative depressions become positive peaks. 



A = and A = ±2 arctan (^^m2/mi^ ; 



A 
A 
A 



±7r and A = ±2 arctan J mi /m2 



±2 arctan A/m2/mi I and A = 0; 



±2 arctan A/mi/m2 ) and A = ivr 



(18) 



For mi = 17, m2 = 83 we have peaks (A, A) = (0,±2.29), (7r,±0.85) and depressions at 
(A, A) = (±2.29, 0), (±0.85, ±vr), where vr — 2.29 = 0.85. These extrema are precisely at the 
positions given by the elements of the density matrix associated with state ( lT3ll . 

The peaks along A = (and ±7r) correspond to phase-diagonal matrix elements which, 
in ( IT6|l . introduce the usual probabilities [1 ± cos A] associated with an interferometer, av- 



eraged oven a random phase A. The extrema centered at A 7^ are phase off-diagonal, and 
directly indicate QIMDS since here the phase state is macroscopic. When A does not vanish, 
probabilities become quasi-probabilities [cos A ± cos A] , which may be negative. 

Now, in ( |T6|) . the rria and m^j dependence is given by the cosine Fourier transform of a 
function obtained by integrating F(A, A) in Eq.( IT711 over A. Because F(A, A) has multiple 
peaks, the final probability contains oscillations as a function of as shown in Fig. H]- if 
we replace the peaks in F{A, A) with 5-functions we recover exactly Eq. f lT4|) . By contrast, 
within SSB the result is equivalent to Eq. (|TU|) . but with A set to zero, which cancels the 
contribution of the off-diagonal peaks. The probabilities then become smooth functions of 
rria, with no dips or peaks; population oscillations disappear. 



FIG. 4: Plot of P{mi,ma) given by Eq. (jl6p versus nia for Nq, = Np = M = 100, mi = 17 and 
m2 = 83. If m2 is even, the central dip is replaced by a peak. 

If we did not count iria and mp but summed over these variables with a given sum 
rria + = N — M, we would get a factor (cos A)'^"'", strongly peaked at A = if - M 
is large. The probability of finding the result set {1711,1712} would then be P(mi, 777.2) ~ 
•I'-TT 1^ + ^'^^ M"^^ ~ '^'^^ A]™^. This still has two peaks in the integrand, which arise since 
the interferometer cannot discriminate between opposite relative phases. But what is now 
obtained is a statistical mixture of these two phases, without any population oscillation; the 
situation is analogous to that described by Eq. (jl]). 

The analysis of Bell violations in Ref. j4] shows that one single missed particle cancels 
the violation. The population oscillations have no special relation to locality, and they are 
more robust. We have shown that, by proper selection of m\ and 7772, one can preserve a 
small central with as many as 5 particles lost. 

In conclusion, two kinds of interference effects occur. One produces the fringes seen in 
the MIT experiments in the merging of two Bose condensates. This effect can be explained 
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by using SSB and either Eq. ([T]) or Eq. (|2]). But the approach using a double Fock state 
preserves a second interference effect: the macroscopic quantum interference that involves 
the off-diagonal elements corresponding to A 7^ 0, and leads to QIMDS that are observable 
via "population oscillations" . 

Laboratoire Kastler Brossel is "UMR 8552 du CNRS, de TENS, et de I'lJniversite Pierre 
et Marie Curie". 
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